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. . . The Haldane-Shastry spin chain can be mapped to the infinite coupling limit of the SU(2) spin 

CO ' Calogero-Sutherland model. We use the Jack polynomials' technology to compute the form 

I factors of the spin operator on the multi-spinon spectrum. The spin structure factor is obtained 

^— «' ■ through a form factor expansion. The expansion is proven to converge in the small momentum limit. 

(N . Numerics based on two- and four-spinons contributions give an approximate result for the infinite 
temperature static and dynamic spin structure factor. 
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I. INTRODUCTION 



Low-dimensional systems constitute fertile breeding grounds for exotic types of physical excitations. Fractional- 
ization of quantum numbers like charge and spin is known to take place respectively in one-dimensional interacting 
^ [ electron liquids and spin chains: in cases such as these, one must forget about weakly coupled particles, and instead 
I ■ adopt a whole new starting point for the description of the strongly coupled physics. Obviously, the identification 

and proper description of this new starting point is often a very involved and risk-prone process. 
^ 1 In this respect, quasi-one-dimensional spin systems have provided one of the sturdiest arenas. Experimental re- 
alizations of systems with fractionalized excitations are numerous and well-documented. Probably the clearest and 
best studied signature comes from neutron scattering experiments on effectively one-dimensional antiferromagnetic 
spin-1/2 chains The excitations seen are not the naively expected spin-1 spin waves, but rather gapless "spinous", 
, which one could losely present as spin-1/2 spin waves. Among many remarkable properties of these excitations are 
^ ' their fractional statistics, intermediate between fermions and bosons, making such a system markedly different from 
O I one obeying convential rules. 

On the theoretical side, strongly-coupled systems like spin chains have in the last few decades presented extreme, 
if not seemingly insurmountable difficulties. The simplest way of explaining this fact might be to say that quantum 
Cn| 1 fiuctuations are very strong in one-d, and cannot be tamed by perturbative approaches. Instead, excitations are 
' strongly nonlinear, and one is faced with the seemingly impossible challenge of either providing an exact solution or 
' risking to miss out completely on the correct physics. 

The quantity of interest to experimentalists (thinking about neutron scattering experiments) is the dynamical spin 
[ structure factor (DSSF). For a chain of N spins at sites Ri, this is defined as 

m: 5"''(q,c.;T) = — -^e^'i-^^^-^') / e-*(5r(0)5f (i))^ (1) 

O . i,3 

; where the angular brackets denote a thermal average. 

a The model of choice for the description of the spin dynamics depends of course on the specifics of the experimental 
I , setup one wishes to describe. The XXZ Heisenberg model often fits the bill remarkably well, at least for very 
'"^ ' low energies. The Bethe Ansatz method could provide most of its thermodynamic properties, but little about 
^ its dynamics. Approximate methods have thus been used to address the computation of the DSSF. The Miiller 
O Ansatz 4j is the best conjecture for the zero temperature structure factor based on exact results and numerics. Finite 
^ ' temperature low-energy features were obtained by Schultz using bosonization 5]. 

^ . To go beyond the field theory limit requires tackling the nonlinear nature of the original model. Considerable 
insight in this direction was provided by the Algebraic Bethe Ansatz method @ and the quantum inverse scattering 
, theory 't'I . Bougourzi et al. used results from the algebraic analysis to compute the exact two-spinon contribution to 
the DSSF of the one-dimensional Heisenberg model 0, 0, ED| • More recently, Maillet et al proved multiple integral 
representations of elementary blocks of the correlation functions j^J. Nevertheless, all these approaches restrict to 
zero temperature and no exact thermodynamic limit is known. The computation of the DSSF at finite temperature 
requires a lot of further efforts. 

However, one of the properties of the Heisenberg model is that the spinous, though dcconfincd, still suffer from 
a residual interaction. The spinous are thus not truly free excitations obeying fractional statistics. There exists on 
the other hand a very convenient alternative approach based on the Haldane-Shastry model, whose Hamiltonian is 

[Hill 



2 



where d{i) = ^ sin It is in the same universahty class as the Heisenberg model, and the long-distance (decaying 
as for large distances) interaction in fact simplifies things considerably: the spinons form an ideal gas of 

particles obeying fractional exclusion statistics The DSSF at T = can in fact be calculated exactly, and is 

given in the thermodynamic limit by 

Sgp Q(^2(g||) - oj)Q{uJ - uJi-{q\i))Q{uj - Ui+{q\\)) 

^g||(7r-g||), wi+(q||) = ^(q|| -7r)(27r-g||), uJ2iq\\) ^ jq\\{27r - q\\). (3) 

This r = formula is made up only of contributions from the two-spinon channel: all higher channels have vanishing 
contributions in the zero-temperature limit. 

One of the very nice features of the Haldane-Shastry model is that its dynamics turn out to be much more easily 
tractable than those of the Heisenberg chain. The Haldane-Shastry model is but the first representative in a wider 
class of solvable models dubbed the SU(N) Haldane-Shastry chains. These are in turn obtainable as a particular limit 
of spin Calogero-Sutherland models, for which an impressive number of exact results are known in the mathematical 
literature. In particular, there exists a Yangian symmetry leading to the identification of a set of eigenvectors 
constructed from Jack polynomials 17]. In short, the whole set of expectation values and transition matrix 
elements one might want to calculate in the Haldane-Shastry model turn out to have a correspondence in terms of 
Jack polynomials. More details on this will be provided in the bulk of the paper. 

Thus, this opens the way to the computation of the DSSF Q at nonzero temperatures for the Haldane-Shastry 
model. Our strategy will be to make use of the technology contained in 17] to compute the form factors involved in 
the spin-spin correlation function needed for the DSSF. This approach has already been used for T — dynamical 
properties of the Haldane-Shastry spin chain [13. [T9| and related models (spin Calogero-Sutherland model [20l] . 
supersymmetric t — J model [2ll l22l]'). The present work is the first to address finite temperature dynamics. 

However for a generic quantum field theory divergences appear when developping a correlation function on the 
Hilbert space. In the context of integrable field theories (ITF), it has been proposed that it could be rewritten as a 
sum free of divergences. The resulting formula, called a 'form factor expansion' (FFE), can be evaluated using the 
scatte ring da ta of the ITF. There is still on ongoing discussion about how precisely the FFE can be implemented in 
ITFs |2a l2ll25ll2ll27Ll2^ . For the case of conformal field theory (CFT), there has been a similar but independent 
proposal for writing finite-temperature correlators in a FFE 29]. It is based on the fractional statistics of the 
quasiparticles building the Hilbert space. Evidence was put forward that it converges quickly to the exact (known) 
result in terms of the number of excited quasi-particles ^29„ .30,] . We will prove in the paper that the latter approach 
applies in the case of the Haldane-Shastry spin chain. 

Our paper is organized as follows. First, we recall all the necessary aspects of the computation of form factors 
for the Haldane-Shastry model using Jack polynomials. We then set out to calculate the form factors themselves, in 
increasing complexity of spinon channels. The form factor expansion is then introduced and proved. We finally put 
the results together to provide an expression for the finite-temperature static and dynamical spin structure factors. 
Discussions and conclusions are amassed at the end. 



^i-(9||) 



II. JACK POLYNOMIAL TECHNOLOGY FOR THE SPIN CALOGERO-SUTHERLAND MODEL 

Using the freezing trick [sj the SU(2) Haldane-Shastry model is obtained by the strong coupling limit of the SU(2) 
spin Calogero-Sutherland model [s^l- The latter describes N particles with coordinates {x^, i — 1 . . . N} moving on a 
circle of length N with the Hamiltonian 

Hs,.ncs- 2l.ax^ + 2A.2|^sin^(^(..-x,)) 

where is the SU(2) exchange between particles i and j. Within the freezing trick, the interaction parameter (3 is 
taken to infinity. The particles are therefore pinned at Xi = i and interact via the SU(2) spin exchange. 

The spin Calogero-Sutherland model is more tractable than the Haldane-Shastry model since it is continuous and 
all his eigenfunctions have been explicitely constructed. We recall in the following the Uglov construction ^3] in 
terms of glj Jack polynomials. The mathematical technology they provide to compute transition matrix elements is 
then introduced. 
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A. Yangian Gelfand-Zetlin basis 

We consider the case N even and N/2 odd, so that the ground state is unique. Uglov determined the so-called 
Yangian Gelfand-Zetlin basis of this model, orthogonal through the Yangian action. They are labelled by the strictly- 
decreasing sequences k = {ki,i = 1...A''}, = {N/2 + 2 — i,i = 1...N} corresponding to the ground state. 
k contains information on both momentum and spin of the excitation. Writing ki — 2ki + ki, ki £ Z represents 
momentum and ki S {1, 2} color^. More precisely, the momentum, the energy and spin of a state described by k are 

i 

Ek^'^-^Y.'^h+P{N + l~2i)/2f (6) 

i 

Sk = \Y.[5k,,2-h^,i] (7) 

i 

For physical applications it is more convenient to work with excitations over the ground state. One identifies the 
state k with the pair 

k = {X {h - kN + i - N, i = I . . .N ^ I), r ^ kN ~ k%) . (8) 

It consists of a zero mode r and a partition A (non- increasing sequence of positive integers) of length iV — 1. The 
Hilbert space is spanned by all possible pairs. We refer the reader to IZj for the expressions of the physical properties 
in terms of (A,r). They will be specified in the next section for the specific case of the Haldane-Shastry spin chain. 



B. Uglov's isomorphism 



Uglov determined an isomorphism Q. between the Yangian Gelfand-Zetlin basis and the Jack Polynomials defined 
through 



fl(fc) = (a;i...x„)'-pf^+^'')({x,J) 



(9) 



where the Jack polynomial P^''^'^ is the limit q = —p, t = —p^, p ^ 1 of the Macdonald polynomial P\{q,t). Then 
one has 



(fc,0(/3,2) -(f^(fc),f^(0)(^,2). 



(10) 



( • J • )(/3,2) is the Yangian scalar product (see jl7i | for a definition) , and ( . , . }/3,2 is the following scalar product in the 
space of symmetric Laurent polynomials 



(/({a^j),5({x.})),,2 = / ^/({^^» 



(11) 



The physical quantity we study in this work is the action of the spin operator. Using Uglov's isomorphism, it is 
given by (we dropped the scalar product indices for convenience) 



(A,r|s±|M,r') = 1 J2^X,r\p2s+i\fi,r')6{Sx ^S^ = ±1), 



(12) 



where Pm = X^i^i ^f' is the power sum symmetric function. One can identify with (s+ + s )/2 and replace the 
6 hy l/2S{\Sx — Sfj_\ = 1). In this paper we only study form factors which satisfy this selction rule. By symmetry 



^ We use here a notation different than Uglov. In fact, this corresponds to the dual representation, called * in his paper. Results are 
not altered, as one has to sum both representations to obtain physical quantities. Our choice for this representation is mostly practical, 
leading to states with positive momentum when ai is positive. 
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between the contributions of pm and p-m (or equivalently the duahty relation (see |17|'). one may finally only consider 
the form factors 



(A,r|s±|M,r')^^ {\\P2s+iWryS{]Sx - = 1). 



(13) 



sez+ 



C. glj Jack technology 



To compute transition matrix elements such as (|13ll . we need some results from the mathematical literature on 
Macdonald polynomials [ssf. 

j 



A, 



FIG. 1: Partition A = (6, 4, 4, 3, 1, 1, 1). 



We consider a generic partition as on Fig. ^ It is a tableau made of cases s labeled by their row number i and 
their column number j. The length (the number of cases) of the row i is A^ and the length of the column j is A^ . One 
then defines 

Definitions for partitions : 

Length: ^A) = max(A^); 
Cardinal: |A| = J2^i '^ote Xh |A|; 

Arm-length: a{s) = A^ — j; 
Leg-length: l{s) = X'j — i; 
Arm-colength: a'{s) = j — 1; 
Leg-colength: Z'(s) = i — 1; 
Content: c(s) = a'{s) — l'{s); 
Hook-length: h{s) = a{s) + l{s) + 1; 
C(A) = {s e X\c{s) = mod 2}; 
H{X) = {s e X\h{s) = mod 2}. 

The quantities of interest in the present work are the norm of a given state, the expansion of power sums on 
Jack polynomials (thus on physical states), and finally the Fieri formula. The latter gives the development on Jack 



polynomials of the product of a Jack polynomial and an elementary function Cr 
available to compute form factors in the space of symmetric polynomials. 
They read 



P. 



(7.2) 



This is the only formula 



Norm 



(7,2)|p(7,2) 



A 



A 



(111) n 



C(A) 



a'{s) + l + j{N -l'{s) - 



n 



a{s) + l+-fl{s) 



1) J-J- a(s) 



l{l{s) + l) 



(14) 
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Expeinsion of power sums 



Fieri formula 



P2s+1= E XAPr''\ (15) 



^ + 7(^(5) + 1) ^ ^^^^ ^ otherwise, (19) 

a(s) + l+7Z(s) ^ " ^ ^ 

with A — /i being a vertical r-strip (at maximum 1 box per row, for a total of r), C^/^ (resp. R\/n) being the 
union of columns (resp. rows) that intersect A — /x. 



III. TRANSITION MATRIX ELEMENTS IN THE HALDANE-SHASTRY MODEL 

We now specify this mathematical background to the case of the Haldane-Shastry spin chain. The /? — > cx) limit 

simplifies a great deal the algebra. The eigenstates described above can be interpreted as multi-spinon states. We 
give their physical properties in the following. Then a general expression for the matrix elements of the spin operator 
is presented along with closed analytical result in the case of few-spinon states. 



A. Spinon interpretation 



In the Haldane-Shatry framework, all physical quantities rewrite in terms of sums over the columns A^- of the 
partition A: 



^ \ - 211] j 



N 



i^-i-y/N)^' 



2w-i ( 2w 



N 



+ (1 - (-)'-)/27V 



[^J for j + r even 
[^] for j + r odd 



(20) 

(21) 
(22) 

(23) 



At the thermodynamic limit (which only is of interest), distinctions between even and odd disappear for the momentum 
and the energy. Defining Xj = irXj/N £ [0,7r], they are 



Pa = nr + ^ Xj 
j 

Ea = Y^Xji-K-Xj) 



(24) 
(25) 



One recognizes the dispersion relation of spinous. An excitation can then be described by a zero mode r and a set of 
particles called spinon defined by each column of a tableau A. We will for the sake of simplicity write such a state 
{rrij = A^ }r. r can be discarded in most of the physical applications. 
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Now we can express all the Jack polynomial technology in terms of the spinons' quantum numbers mj. We will use 
the short cut notations 



r(Lm/2j+l)r(l/2) ^ Nx 

^-"^('"^ ^ r(LW2j + i/2) " V — ^'^^ 

7odd(m) = - ^-J (27) 

27even(TO) TT V 2 



The norm is 



^ ^ {P{m,}\P{m,}) _ -A- 7»+l(^) TT 7»-j(™» ""^j) _ TT / 27r 

One shows that the power sum operators decompose into 2-spinon states (mi, 7712) such that mi + TO2 = 2s + 1 
with 

X(mi,m2) = (-)'"'7o(mi - m2). (29) 
Specifying the Fieri formula to /i = (mi, . . . , m„), A = (pi, . . . ,Pn+i) — with possibly Pn+i =0 — it gives 

^/ ^ TT l^-jiP^-Pj+lh^-3i"^'^ ' ™i) /gg) 



B. Matrix elements 



Different strategies apply to the evaluation of the transition matrix elements (A|p2s+i|m)- The first one is to write 
/i and A into elementary functions e by inversion of the Fieri formula p7(l , use the decomposition p5(l , and apply the 
Fieri formula successively on the multi-e state. This solution is quite unpractical, because inverting the Fieri formula 
becomes increasingly difficult with the number of spinons considered. It is trivial for 1 spinon, and leads to a rather 
cumbersome expression already for 2 spinons. Still, we can use it to rewrite the power sums as a sum of 2-e states. 
The result is interestingly simple 

P2S+1 = I i2^-r^^^' (2r + l)e,+,+ie,_, (31) 

r=0 

To prove it, one uses p7|l on H31|l . it leads to (|29|l thanks to the equality 

It gives way to a generic strategy: use decomposition ()r>l(l and use the Fieri formula twice. The result is (with 
normalized states) 

(Ab2.+i|^) = I + 1^1 - 2W\MvK\, \^ (33) 

V 

with |A| - 2s + 1. 

We will now evaluate the contribution of several channels. Eq. doesn't give a closed analytic expression for 
the form factors. Only in a few cases can it be so reduced. The results may be conjectural, then confirmed through 
the following sum rule 



"S^ {^J■\pmP-,n\^J■) — — >\p.\ (34) 

which can be proved by simple Jack polynomials' algebra. 
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2 spinons This is the only channel present at zero temperature. It has been conjectured by Haldane, then proved 
using the simplectic ensemble P^], and finally Yamamato et al. |0| obtained it using Uglov's technology. 

We recall the result and give the thermodynamic limit 



u M / 7o("ii - m2)7i(mi - m2)7o[iV]7i[iV] / 7r(a;i - X2) 

y 7i(mi)7o(A' - mi)7o(m2)7i(iV - TO2) V ^JE{ml)E{m2) 



1 spinon 



((m )\pm'-„i\{m)) = W— 7T7 ^ ^ —( ^ (36) 

y 7o(to)7o(A' - m') \x[t:-x')J 



2 spinons 

((TOi,m2)|pm|(mi,m2)) 



' 70(771^)71 (m^i)7o(A^ - TOi)7i(A^ ~ m2) 
70(^2)71 (mi)7o(A^ - TOi)7i(A^ - m^) 



/ 7o(fc)7i(fc') -r I _ 

7o(fc')7l(fe) ItTOi-mi, 



7o(fc')7i(fc) :f™/ _ ™„ ("^7^ 



7o(fc)7l(fc') "'2 

' Zt)lt') G{m,l = m'2-m2,k,k') otherwise. 

2 I ^ I 2 I ^ ^ L 2 -I ' 2 '' V I 2 I ' 2 



G(m, fc, fc') = y (_)»('=+'+i) IiW ill 



^^1) r([i]) r(L^j + i) r([|] + i) 



7o(*) r(L^j + i) r(L^j + i) r(L|J + i) r(r^l , 



(38) 



with m — m'^ + m, — mi — m2, fc = mi — m2, k' — m'l — mj. This can be interpreted as an SU(2) generalization 
of result (28) of for g = 1/2 quasiparticles. 

Higher channels It is generally not possible to obtain a closed analytical expression for the other channels, except 
in a few cases. Such are, for example, [(mi) — > (mi,n, n')] and [(mi) — > (71, n',mi)] which equal [0 — *■ (n, n')]. 
Nevertheless, formula (|33|l can be used to give exact numerical results. Getting the thermodynamic limit directly 
is a challenging but fruitful work. We leave it as an open question. 



IV. CORRELATION FUNCTIONS AT FINITE TEMPERATURE 



Before addressing the computation of the DSSF, general considerations on finite-temperature correlation functions 
are needed. For a local operator O, it is 

i:„EA,.(A„|0|A„)exp(-/?£;,,.) 
E„ Ea„ exp(-/3SA J 

where n is the number of quasiparticles, and A„ a state with n quasiparticles. 

As recalled in the introduction, divergences appear in the correlators of the right-hand side that need to be re- 
summed. For ITFs, LeClair and Mussardo proposed such a resummation as a form factor expansion on the basis of 
the asymptotic particle states in the zero-temperature theory |^ 

n 

{0)t ^Y.J2(^^\'^\^^)"-]1''t{ExJ (40) 

n Ati i—1 

The irreducible form factor (A„|0|A„)irr is obtained thanks to the Form Factor Bootstrap (FFB), and ut^E) is the 
filling factor determined by the Thermodynamic Bcthe Ansatz (TBA). 

In the following, we show that the FFE apply also for the Haldane-Shastry spin chain in the thermodynamic limit. 
First we obtain the thermodynamic properties of the spinons, then give the expression of the irreducible form factor. 
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A. Exclusion statistics 

A form factor expansion similar to H40fl was proposed for CFTs . It led to identify a two-body S-matrix for the 
CFT in the thermodynamic limit 

S = exp[2iTr{6 - K)@{0)] (41) 

where K is the exclusion statistics' matrix of the quasi-particles of the theory (which play the role of asymptotic 
states) . 

Fractional exclusion statistics is a tool introduced by Haldane il5||| for the analysis of strongly correlated many-body 
systems. It is only based on the assumption that the Hilbert space is finite-dimensionnal and extensive, i.e. particles 
are excitations of the considered condensed matter system, so it is a very generic concept. The statistics are encoded 
in a matrix K ~ [Kij ) correponding to the reduction of the available Hilbert space for particle of type i by filling a 
one-particle state by a particle of type j. This is then a generalization of the Pauli principle. 

For spin- 1/2 spinons with species i — ±, the statistical matrix is [s^ 

K=(^ \ I ) (42) 

As such, their 1-particule distribution functions generalize the familiar Fermi-Dirac and Bose-Einstein ones. They are 
derived from 1-particle grand canonical partition functions Gi given by the lOW equations [35„ ,36,. „37j 

^)X{Gf--z. (43) 

where Gi depends on the generalized fugacities Zj — e''^^^^^^. The one-particule distribution functions are obtained 
through 

n,{e) = z~\og\[G, (44) 

' 3 



In our case, where species are i — (+, — ), we obtain 



from which one deduces that at zero magnetic field, these distributions match the Fcrmi-Dirac distribution function. 
It means that spinons can be considered as having fermionic statistics. This is indeed what is done when they are 
labeled by ordered numbers within partitions, the spin degree of freedom is hidden. The norm formula (|28|l shows 
that for a multi-spinon state, the labels have to be strictly decreasing, as for spinless fermions. 



B. Irreducible form factors 



For ITFs, the definition of the irreducible form factors comes a priori from the FFB. Within our framework (based 
on pgf ') it is not even necessary. Their definition is 

(A„|0|A„) = (A„|0|A„)irr + 5I(^I^^K>"T (46) 

with An being a sub-state of An (a substate being a state where some of the spinons have been taken out). Nonetheless, 
the FFB insures that irreducible form factors don't carry divergences, which we can't prove here. This calls for further 
understanding of form factors in fractional statistics' theories. 

Here follows a sketch of the proof that the irreducible form factors (|46|) give the correct FFE 14U|) . It only uses 
the fact that the colorless spinons are free fermions, which shows up by the strict ordering of their quantum numbers 
within a multi-spinon state. As such, the proof is similar to |38j |. but in the discrete case. 
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We first remark that (OlOlO)!!!' trivially comes with the factor 1. Let us consider the factor of £)(m) = {{m)\0\{m)) 
The contribution C„ from n quasi-particles is obtained recursively, isolating the spinous whose label match with m: 



n 

^" = 7 E E ^(™«) n exp(-/?i?™. ) 



rni>...>?nn i 



1=1 



ri-1 



E n --'"-^ E ^(-)e-^^'"(l - E '^'"-.) 

7ni>...>m„_i i— 1 m -i— 1 

n-1 



(47) 
(48) 



E n^^''"-E^( 

mi>...>mii_i i=l m 
n-2 



me 



n-2 



E n 

m\> ...>mn-2 i— 1 



5]Z?(m)e-2/'^'"(l-^5„,„ 



i 

^ []exp(-/3i?™J. 

mi>...>mi j — 1 



Then summing over n gives 



= Y,D{m)e 



-PE, 



,-iPE„ 



1=0 - 



i=Q 
1 



exp (/?£;„) + 1 



(49) 

(50) 
(51) 

(52) 
(53) 
(54) 
(55) 



The demonstration follows the same lines for the higher form factors. Eq. 0U| is quite convenient. Indeed, the 
denominator is suppressed, so it can be seen as a perturbative series, provided that it's converging. It also has a clear 
physical interpretation with quasi-particles and filling-factors. Now it can be applied to the particular case of the 
dynamical spin structure factor (DSSF). 

V. FINITE-TEMPERATURE SPIN STRUCTURE FACTOR 

Now having all the methodological ingredients we address the main subjet of our paper: the computation of the 
finite-temperature spin structure factor of the Haldane-Shastry spin chain. We first recall the zero-temperature result. 
It was first obtained by Haldane and Zirnbauer |16| using the supermatrix method. It can also be accessed with 
Jack polynomials flsL ll^ . Only the [0 2] channel contributes so it is straightforwardly 



So(q,uj) — / dxi , — ^ , 



dx2 . „, , '^(^1 +X2 — q)^ (Eixi) + E{X2) — ijj) 



(56) 



At finite temperature, all channels contribute, one needs to sum them all through a FFE. We detail in the following 
the static structure factor S't(<?) — J dio ST{q,^^) and the DSSF. 



Static structure factor 



The static spin-spin correlator can be described by a one-point function 

Sxiq = 2tts/N) = {-^P-{2s+i)P{2s+i))t 



(57) 
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(this is not the exact one, but equal at the thermodynamic hmit) To compute it, we use the 1-point FFE and the 
transition matrix elements obtained before. 

Due to the complexity of the form factor and the lack of FFB, it is hardly impossible to obtain analytical results 
on the correlation functions. The only point for which we could obtain the correlator exactly is q = 0. We compute 
it as the thermodynamic limit {N — > oo) of (|57l) at s = 0. The FFE is 



{p~iPi)t = ^^nF{E^){fi\p^ipi\^i)i„ (58) 

(59) 

n 

mi, . . . ,mn}\p-ipi\{mi, . . . ,m„})irr (60) 



N 

oo 
n=l 



In the appendix we show that 



from which wc conclude 



= / dx[nF{E{x))r (61) 







dxe"''^(") (63) 



This result brings strong evidence of the power of the FFE to obtain finite-temperature correlation functions. This is 
the main achievement of the paper. Let us remark here that Fn is the contribution of the n-spinon states in the FFE. 
They appear to be of the same order in temperature, meaning that the FFE is not a low-temperature expansion^. 
It really is a perturbative expansion, in the sense that the different contributions decrease exponentially with the 
number of spinous involved. 

Nonetheless, St{<1 = 0) is a low-energy feature of the theory that can be obtained with a bosonization approach. 
As in '3^, we observe that the FFE is mostly powerful in regimes where other simpler methods apply. To calculate 
the whole static structure factor, we then rely on numerics. We work in the infinite-temperature regime to compare 
the FFE with the expected result. In this limit, correlations are only local and we expect the static structure factor 
not to depend on the momentum and have the value St=oo(<1 = 0) = 1. 

To observe the FFE perturbative power, we studied the contribution from 2- and (2+4)-multispinon states. What 
we understand here as 2-spinons is the sum of the contributions of the [0 <-> 2] and [1 ^ 1] channels, and as 4-spinons 
the [1 ^ 3] and [2 ^ 2] channels. Results are gathered on figure [3 At first we observe the convergence to be rather 
good. This strongly supports our approach. Still problems arise in the vicinity oi q = tt. Finite-size effects starts 
appearing at the level of 4-spinons contributions, and increase with a higher number of spinous. The only way to 
tackle them is to have an enormous computing time. We admit this is the most important weakness of the approach. 
This emphasizes the need for a direct thermodynamic limit computation. 



B. Dynamic structure factor 



We show on figure |31 the sum of the contributions from 2- and 4-spinons to the infinite-temperature DSSF. As in 
the previous paragraph, finite-size effetcs clearly appear around q = n. This makes the comparison to experiments 
hardly possible. We still note that some basic features, like the double arch shape, agree qualitatively. This feature, 
already present at zero temperature, comes from form factors where the spin operator acts on only one of the spinon 
in a multi-spinon excitation, leading to the spinon dispersion relation. We leave more comments to the concluding 
section. 



^ This would have been expected on the ground that more thermal energy is needed to excite more particles. This is not true anymore 
for massless particles. 
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FIG. 2: Static spin structure factor. 2 is the contribution from 2-particle form factors with A'^ = 10000; 2+4 is the contribution 
from 2- and 4-particle form factors with A'^ = 500; the two different curves are obtained with two labeUings of the spinous: for 
the upper curve, < m < N (the normal one), for the lower curve, ^ m ^ N; they should match in the A'^ —> oo-limit (up to 
a 5-peak on g = vr), so we expect the thermodynamic limit to lie inbetween the two curves. 

VI. CONCLUSION 

In this paper we studied the finite-temperature spin-spin correlation function of the Haldane-Shastry spin chain. 
We used the multi-spinon basis obtained as the infinite couphng limit of the Yangian basis of the spin Calogero- 
Sutherland model. Form factors of the spin operators were computed thanks to the glj Jack polynomial technology. 
These form factors were gathered within a form factor expansion of finite-temperature correlation functions to give 
physical quantities directly comparable with the experiments. Though the Haldane-Shastry is not a quantitatively 
appropriate model to describe real spin chains such as KCuFa, it should at the qualitative level. As any finite- 
temperature properties for the Heisenberg model are yet beyond reach, serious theoretical insights are brought by the 
Haldane-Shastry model. We recall that it is simpler due to the 1/r^ interaction responsible for the thermodynamic 
freedom of the spinous. 

Our main analytical results are: the formal expression of any form factor of the theory, a closed analytical expression 
for the simplest of them, the proof of the finite-temperature FEE, and the value at zero-momentum of the static 
structure factor. Along with numerical results, they make evidence of the EFE technique to approximate finite- 
temperature correlations. As a drawback, it demands huge computation power. This calls for further theoretical 
refinement. 

Work is on progress to provide the Haldane-Shastry model with integrability features such as the FEB. Within 
such a description, thermodynamic irreducible form factors would be obtained much easier using scattering proper- 
ties. Problems to develop this approach come from the intrinsic discrete nature of the Haldane-Shastry spinous. But 
necessary effort has to be made to gather comprehensively ITEs and CETs in a common framework. Further under- 
standin g is also needed in the link between the Haldane-Shastry model and the symplectic random matrix ensemble 
used in |lg. The analysis performed at zero temperature could be extented at finite temperature. It would also be 
interesting to treat the finite-t emp erature dynamics of other inverse-square interaction models and compare it to their 
zero-temperature exact result |2Cl |22|. It comes as a simple generalization of the method used in the present paper. 

At the numerical level, finite size calculations can be performed. For sufficiently small sizes, the exact correlation 
functions are accessible using the analytical results of this paper. The FEE can not be used, for it is based on the 
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FIG. 3: Dynamic spin structure factor; the contributions from 2- and 4-particle form factors with A'^ — 500, in absolute value. 

thermodynamic limit, rather direct apphcation of l|39|) is necessary. Using directly the Yangian multiplets is another 
solution, as in 39] for zero temperature. 
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APPENDIX: ZERO-MOMENTUM LIMIT OF THE STATIC STRUCTURE FACTOR 

In this section we show that 

F„ = ^ ]^n_F(i;,„J({mi, . . . ,m„}|p_ipi|{mi, . . . ,m„})i,.r = / dx [np{E{x))]" (A.l) 

It is clear from the conservation laws that the intermediate states for each form factor is only different from the initial 
state by an increase of 1 of one of the spinons' m. Quite generally, we thus need 



|({toi, . . . , mfc + 1, . . . ,to„}|pi|toi, . . . , mfc, . . . , m„})|^ ~ Y[ ot^^y^kHriik - rui) 



a_(_(— m) 



"foim)"fi{Tn + i; 
Jiim)jo{m + 1] 

7i(m)7o(m - 1) 
7o(m)7i(m - 1) 



(A.2) 

(A.3) 
(A.4) 
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Explicitly developing the irreducible form factors and separating different affected m's, we obtain 

n 

Fn= Yl U.^f{e^.)J2f!: (a.5) 

mi>--->m„ i=l k 

= E n (^^ - mi) (A.6) 

JC{l,...,k-l,k+l,...,n} iGJ 

J>C{fc+l,...,n} ieJ> J<C{l,...,fc-l} iGJ< 

with dj{k,i) the distance between /c and i in the subset J. Thus we can write = F^^F^^. 

We will know obtain J2{mi>k} ™ ^ recursive way (the proof for follows the same lines), putting k = 1 
without loss of generality. We first perform the sum over m„. Writing 7 = {2,...,n— 1}, one separates the ensemble 
of J as 

{ J C / U {n}} = { J + J U {n}, J Cl} = {J + JU {n}, J Gl, \J\ even} U { J + J U {n}, J C J, | J| odd}. 
In the first (resp. second) subset, the distance in J between 1 and n is odd (resp. even). This proves that 

F^> = {at{m, - m„) - l)F^>r^ + (a+(mi - m„) - l)F^>f'' (A.8) 

where the superscripts even and odd corresponds to restrictions of the expression of -F^^i to J subsets with even or 
odd cardinal. Repeating this recursion one ends up with 

Fn^ = F^>r"[(a+(mi - m„) - l)at (mi - m„_i) + 1 - at(mi - m„)] 

+ F^>f'^[{a±{mi - mn) - l)a+(mi - m„_i) + 1 - a+(mi - m„)] (A.9) 

Then one easily shows that 

[(a'|(mi - m„) - l)at{nii - m„_i) + 1 - atimi - mn)] ^ at (mi - m„)(5„i„_„i„_i (A. 10) 

[(at (mi - m„) - l)a|(mi - m„_i) + 1 - a|(mi - m„)] ~ a|(mi - mri,)Sm„-m„-i (A-H) 

thus proving the recursion 

F}^^ = F^>,6m„-m^., (A.12) 
As this is true only for n > 1, we obtain for the following 

F^ = {at{mi - ms) - l)Sm2,...,m„ (A.13) 

pkyil,n ^ (Q,+ (TOfe - mfe+i) - l)(a:,:(mfe - mfe_i) - l)^mi,...,mfe_i^mfe+i,...,m„ (A.14) 

F^ = (Q!+(m„m„_i) - l)6rm,...,mn-i (A.15) 
Now, one has a't{mk — TUk+i) — 1 = Smi,,mk+i ^^'^ Q!^(mfc — mfe_i) — 1 = so that 

F^ = 6k,l6m„...,m„ (A.16) 

Finally 

N-l 

= E lME{m))r = / dx [nF{E{x))]^ (A.17) 

m=l "^0 

which ends our proof. 
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